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We analyze the evolution of cosmological perturbations in the cyclic model, paying particular 
attention to their behavior and interplay over multiple cycles. Our key results are: (1) galaxies 
and large scale structure present in one cycle are generated by the quantum fluctuations in the 



■ preceding cycle without interference from perturbations or structure generated in earlier cycles and 

, without interfering with structure generated in later cycles; (2) the ekpyrotic phase, an epoch of 

gentle contraction with equation of state w 3> 1 preceding the hot big bang, makes the universe 



homogeneous, isotropic and flat within any given observer's horizon; and, (3) although the uni- 
verse is uniform within each observer's horizon, the global structure of the cyclic universe is more 



complex, owing to the effects of superhorizon length perturbations, and cannot be described in a 
uniform Friedmann-Robertson- Walker picture. In particular, we show that the ekpyrotic phase is 
so effective in smoothing, flattening and isotropizing the universe within the horizon that this phase 

X. 

alone suffices to solve the horizon and flatness problems even without an extended period of dark 
energy domination (a kind of low energy inflation). Instead, the cyclic model rests on a genuinely 
novel, non-inflationary mechanism (ekpyrotic contraction) for resolving the classic cosmological 
conundrums. 
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I. INTRODUCTION 



Cosmological observations support the idea that the part of the universe we now observe 
emerged from a hot, radiation dominated and expanding state. The original hot big bang 
theory did not specify how this radiation was generated. Inflationary theory [l| postulates 
a phase of superluminal expansion, driven by scalar field potential energy which ultimately 
decays into radiation. By contrast, in the ekpyrotic [2j model, the radiation is generated 
by a brane collision, following an earlier empty phase. The earlier phase is contracting 
from the viewpoint of Einstein-frame four dimensional effective theory. A transition from 
Einstein-frame contraction, to expansion is also invoked J,n_the pre-big bang model jjj . 



Both the inflationary 



plausibly give a spectrum of 



and ekpyrotic models 

primordial fluctuations that is nearly scale invariant and adiabatic, in accordance with recent 
cosmic microwave background 9] and large scale structure jl(| observations. Indeed, these 
two theories seem to be the only ones that are generically able to give such predictions 11 1. 
However, both theories are incomplete. For inflation, it seems inevitable that an initial 



singularity is still present 



12j (but see 13|,llJ]) only a finite affine parameter distance away, 



meaning that initial conditions at the singularity are potentially significant. In the ekpyrotic 
scenario, the initial emptiness of the contracting phase is not explained. 

The cyclic model \ioL ll(| builds on the ekpyrotic one, essentially stacking a whole series 
of ekpyrotic histories together. One post-bounce expanding phase links onto another con- 
tracting phase. This contracting phase passes via a bounce into a new expanding phase, 
and so on. The universe evolves cyclically, so that the puzzle of a "fundamental beginning 
of time" is at least deferred into the very distant past, and perhaps avoided altogether. 

The ekpyrotic and cyclic models both have a higher dimensional interpretation, inspired 
by M-theory. This perspective is critical for understanding how perturbations can propagate 
through a bounce, and much recent work has focussed on this issue |8, 




four dimensional perspectives on the matching see 
Here we address other critical issues for perturbations in the cyclic model, and are able 
to use a four dimensional effective description almost throughout. We simply assume the 
essential features of the five dimensional matching prescription are correct, and apply them 
to the four dimensional effective theory. 

Our analysis shows that the galaxies and large scale structure in any given cycle can 
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be generated by the quantum fluctuations in the preceding cycle without interference from 
perturbations or structure generated in earlier cycles and without interfering with structure 
generated in later cycles. The global structure of the cyclic universe is more complex: al- 
though the universe can be described as a nearly homogeneous and isotropic within any 
observer's horizon, the global structure cannot be characterized by a uniform Friedmann- 
Robert son- Walker picture. Our results further show that the ekpyrotic phase alone is suf- 
ficient for resolving the horizon and flatness problems and that an extended phase of dark 
energy domination or any other form of inflation is completely unnecessary. This makes it 
clear that the cyclic model is a genuinely novel, non-inflationary approach to cosmology. 



II. OVERVIEW OF THE CYCLIC MODEL AND ITS PERTURBATIONS 

The cyclic model assumes that we live on a brane in a special configuration of a higher 
dimensional theory such as M-theory. Away from a bounce, the universe can be treated 
using a four dimensional effective theory consisting of gravity coupled to one or more scalar 
fields. Assuming the background universe is spatially flat, the metric is 

ds 2 = -dt 2 + a 2 (t) 5 ij dx i dx j (1) 

where a(t) is the scale factor. The main imprint of the higher dimensional theory on the 
effective picture is through the addition of one or more scalar fields (f> with a potential V(4>). 
This potential performs many functions in the cyclic model, including that of describing the 
dark energy responsible for the cosmic acceleration observed today. Through most of this 
paper, we shall describe the model in terms of a single scalar field, although we note that 
generically more than one scalar field is involved. The scalar field <fr satisfies 

4> + 3H<P = -V >4> , (2) 

in the background (Q), where dots denote derivatives with respect to t and H = a /a. 
Ignoring, for simplicity, the coupling between ordinary matter and 0, the Friedmann equation 
is 

H 2 = \(p+\tf + V{^ (3) 

in reduced Planck units (87rG = 1), where p is the energy density of ordinary matter and 
radiation. 
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FIG. 1: An example potential V((p). This plot shows where <j) is on its potential at each stage in 
a cycle. The equation of state parameter of the background solution is denoted by w. 

The potential V(4>) is chosen by hand at present, but should ultimately be derivable from 
the higher dimensional theory. It must be of a certain form in order for the cyclic model to 
work 28j . A useful potential with the desired properties is 

V^) = V {e b *-e-^)F(4>) (4) 

(see Fig.HJ). Here Vb is of order today's dark energy density, b is non- negative (and typically 
<C 1) and c is positive (and typically ^> 1). F(<ft) is a function whose precise form is 
unimportant, but which tends to unity for greater than end and to zero for less than 
</> en d. The resulting potential V((f>) has a large negative minimum, denoted V en d, at (/'end- 
While the explicit exponential form used here is convenient for analysis, note that the cyclic 
model actually works for a very wide range of potential forms, the only conditions being that 
they have a steep, negative and strongly negatively curved region over the observationally 
relevant range of the scalar field. 

Of central importance to the cyclic model is the ekpyrotic phase, in which the universe is 



slowly contracting and the scalar field is rolling slowly down its steeply declining, negative 
potential. For our example potential, the negative exponential dominates, V(<p) ~ —Voe~ c ^, 
and the background universe enters an attractor scaling solution, 

a(t) oc (-t) 2/c2 oc e^ /c , H = oc -e"^ 2 , w « c 2 /3 > 1, (5) 

in which t is negative and increasing, and w is the ratio of the pressure to the energy density. 
Notice that, since c ^> 1, as the scalar field moves over a substantial range in Planck units 
towards en d, the scale factor contracts by only a modest factor. In contrast, the Hubble 
parameter H grows dramatically, beginning from values comparable to today's value, and 
growing to values corresponding to high energy scales. 

The scaling solution is only relevant as long as > en d, and the function F(<f>) is 
effectively unity. Once passes the potential minimum, the potential energy is quickly 
converted to kinetic energy and the solution enters a kinetic energy dominated phase, with 

a(t)oc(-*)5<xe* /V5 , H = ^- oc -e^v 7 ^, w ~ i. ( 6 ) 

Ob 

When lifted to higher dimensions, this solution describes two colliding branes (one with 
positive tension and the other with negative tension), whose scale factors remain finite 
even as the four dimensional scale factor a(t) tends to zero and the scalar field tends to 
— oo. Near the collision, the four-dimensional Einstein-frame metric and scalar field become 
singular coordinates; however, five- dimensional quantities like the metric on each brane, 
and the inter-brane distance, are perfectly finite. The matching of perturbations across the 
bounce is therefore performed within the higher dimensional setting. 

As the branes emerge from the collision, the solution followed is nearly the exact time- 
reverse of (jHJ); the radiation and matter produced at the bang and a modest enhancement 
of the kinetic energy of have a negligible effect while < en d- There is a brief w ^> 
1 expanding phase after passes en d moving to positive values, but the excess kinetic 
energy in quickly overwhelms the potential energy V(<ft) and the universe enters a second 
expanding kinetic phase (Figure 1). As is shown in the Appendix, the expanding w ^> 1 
phase is of modest duration and for the remainder of this paper it can be safely ignored. 
It is then convenient to describe all three kinetic energy dominated phases, namely the 
contracting and expanding kinetic phases with < en d, and the second expanding kinetic 
phase with > en d, as a single kinetic phase, and we shall generally adopt this terminology 
throughout the remainder of this paper. 
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As we continue into the expanding phase, the kinetic energy in redshifts away as a~ 6 
and the universe becomes dominated by the radiation that was produced at the bounce. The 
net expansion in the entire kinetic phase is ~ e 27//3 , where 7 = ln((— Knd) 1 ^ 4 /^), and T rh 
is by definition the temperature of the radiation when it comes to dominate. As shown in 



Ref. 28], cyclic models require 7 ~ 10 — 20 in order to be compatible with observation. The 
additional Hubble damping due to the radiation has the effect of slowing down to a halt 
on the positive potential plateau. Then, the scalar field begins to gently roll downhill. The 
matter era passes and the universe enters the dark energy phase. Eventually, the rolling of 
carries it off the plateau. The accelerated expansion due to dark energy reverses to slow 
ekpyrotic contraction. The universe heads towards the next bounce and the next cosmic 
cycle. 

The evolution of the background universe is illustrated by the inner, solid, track of the 
"wheel" of Fig. |21 Fig. ^ has also been labelled to show where is on its potential at each 
stage in the cycle. See also Fig. 0] for a summary of the behavior of key quantities. 

Note that while the bounce itself is nearly symmetrical, the background evolution for 
> 0end is highly asymmetrical, and the scale factor undergoes a large net expansion from 
cycle to cycle. As explained above, the kinetic phase gives a net expansion of |7 e-folds. The 
ensuing radiation phase gives a large number of e-folds of expansion, and the matter phase 
adds a few more. We may approximate the combined number from the latter two phases as 
iV ra d = ln(T r h/T ), where T is the cosmic microwave background temperature today. Dark 
energy adds another potentially large number of e-folds iVdark- By contrast, in the ekpyrotic 
contraction phase, the scale factor contracts by a very modest factor (from Eq.EJ a oc 
So there is a large net expansion every cycle of approximately 27/3 + A^ rad + iV dark e-folds, 
which is critical for the fate of the model when perturbations are considered (see Sec. IVI|) . 
The large net expansion also plays a key role in diluting the entropy density from cycle to 
cycle, and, as we shall see in Sec. IVHH in the cyclic model's solution to the flatness, isotropy 
and horizon puzzles. 

While the scale factor grows with each new cycle, locally measurable quantities like 
the Hubble parameter and the density undergo periodic evolution. The Hubble parameter 
decreases by 27 e-folds in the second kinetic energy dominated expanding phase, and by 
2iV rac i e-folds in the ensuring radiation phase. By contrast, in the ekpyrotic contracting 
phase, H increases in magnitude very rapidly, by a total of iV e kp = In (^\/—V cn d/Voj e-folds. 



Match 




FIG. 2: A "wheel" diagram indicating the behavior of both the background solution (inner, solid, 
line) and perturbations (outer, dashed and dotted, lines) in the cyclic model. "Start" marks the 
point in the cycle at which our perturbation analysis in Sec. 11111 begins. 

Since V ~ Tq in order of magnitude, we find Nekp ~ 2(iV rac i ~\~ / y'), which is also the condition 
that the Hubble constant returns to its original magnitude after a cycle. 

In this paper, for the study of perturbations, we shall consider a universe that contains, in 
addition to the scalar field <f), cold dark matter and radiation, which we treat as perfect flu- 
ids. This description captures the broad features of cosmology well enough for our purposes. 
The metric has scalar, vector and tensor perturbations. We assume that the perturbations 
can be well treated in linear theory, in which case the three sectors decouple and we can 
follow perturbations Fourier mode by Fourier mode. We ignore the vectors and tensors, and 
concentrate on the scalar sector, in which the matter density perturbations exist. We work 
in longitudinal gauge, in which there are no time-space or traceless space-space perturba- 
tions. Perfect fluids and scalar fields do not support anisotropic stress, so the gravitational 
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potentials are equal, and the perturbed metric reduces to: 

ds 2 = -(1 + 2®)dt 2 + a 2 (t) (1 - 2^)5 ij dx i dx j J (7) 

where $ is the Newtonian potential. 

Within any one cycle, the Einstein and matter equations fully determine the classical 
evolution of the perturbations. This classical evolution is indicated by the outer, dashed 
loop of our wheel diagram (Fig. EJ. 

The cyclic model, like inflation, relies on the amplification of quantum fluctuations to ini- 
tiate structure formation. It is helpful to think in the Heisenberg picture. Here the quantum 
field mode operators satisfy the classical equations of motion but even if the quantum ex- 
pectation value of a mode amplitude is zero, the expected variance cannot also be zero (just 
as for the ground state of a simple harmonic oscillator, for example). When the equation 
governing a mode moves from having oscillatory solutions to growing and decaying solutions, 
the quantum variance will also grow, as the square of the classical growing mode amplitude. 
As far as the evaluation of future expectation values is concerned, it now becomes possi- 
ble to accurately approximate the quantum picture with a classical one in which the mode 
amplitude is treated as a random variable with mean and variance given by the quantum 
calculation. We say that a perturbation has been generated when the classical probabilistic 
description becomes accurate. In the cyclic model perturbations can be generated during 
both the dark energy and ekpyrotic phases, and quantum fluctuations in one cycle become 
classical stochastic perturbations by the next. Of course, this stochastic contribution to a 
mode amplitude is only important if it is comparable to or greater than that which is already 
there from the classical evolution. "Quantum generation" of perturbations is represented by 
the dotted line in the wheel diagram, Fig. |21 

To complete the perturbation loop in our wheel diagram we need to know how to match 
perturbations across a bounce. This is the one place where the four-dimensional effective 
picture becomes invalid and results obtained in higher dimensions must be used. Recent 



work 
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29] suggests how this occurs, with long-wavelength growing modes going in to 
the crunch matching onto growing modes going out from the bang. This matching occurs 
in a manner that, for long wavelengths, is independent of wavelength. 

Hence we are now able to follow perturbations through multiple cycles in the cyclic model, 
allowing for both quantum generation and matching across the bounce in addition to the 
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classical evolution. 

III. QUANTUM GENERATION OF PERTURBATIONS 

In principle, when describing a cyclic model, one can start anywhere on the wheel diagram, 
Fig. El For simplicity, we start well into a long-lasting dark energy phase in which the 
universe has become very homogeneous and flat and any pre-existing matter, radiation, or 
scalar field perturbations have been redshifted away to negligible levels. Later on, in Sec. lVll 
we show that the consistency of the model does not require that the dark energy phase be 
long-lasting. 

In a universe containing only a scalar field 0, the Einstein equations fix the scalar field 
fluctuation 5<fi in terms of the Newtonian potential $ and its time derivative: 



Thus there is only one true scalar degree of freedom, which we take to be <3>, and $ satisfies 
the second-order differential equation: 



This equation can be used both for the classical evolution of $ and, as discussed above, 
for determining the variance of the fluctuations generated quantum mechanically. Since, by 
assumption, there are initially no classical perturbations, we turn to the quantum case. 

For the quantum calculation we need to pick a suitable quantum state for each Fourier 
mode. Just as in inflation, we assume that when the evolution of a perturbation mode is 
"gradient dominated" (sometimes called "subhorizon" ) , the scalar field fluctuation 5<f) is in 
the appropriate incoming adiabatic vacuum state. Eq. (JBJ) is then used to determine the state 
of $ in this period (see e.g. J?]). We then evolve forward in time using (0) until the spatial 
gradients become negligible in the time evolution. Now, the mode is said to evolve in an 
"ultralocal" (sometimes called "superhorizon" ) manner. The modulus squared of the mode 
amplitude then gives the quantum variance which, when the quantum picture is replaced by 
the stochastic classical one, becomes the classical power spectrum on that scale. Repeating 
the calculation for different comoving wavenumbers allows us to build up the complete power 
spectrum. The power on a given scale changes with time in accordance with Eq. Q, but all 
modes that are in the long- wavelength, ultralocal regime will evolve in concert. 



ZSff, = $ + H$. 



(8) 
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We perform the above procedure, solving Eq. Q mode by mode with the appropriate 
initial conditions, to build up a power spectrum for the perturbations. All modes of interest 
start off gradient dominated and end up in the ultralocal regime. Longer wavelength modes 
start to follow ultralocal evolution sooner, shorter wavelength modes later. The very shortest 
wavelength modes go ultralocal only in the kinetic energy dominated phase just before the 
big crunch. Once all modes have gone ultralocal, the whole power spectrum evolves in 
concert and simply grows in amplitude as the bounce is approached. 

The detailed shape of the power spectrum depends on the exact background evolution 
and the specific details of the scalar field potential V ((/)). Fig. El shows a power spectrum 
for a typical model (with b = 0.1 and c = 30), evaluated close to the crunch when all the 
modes are evolving ultralocally. The comoving wavenumber is denoted by k. Observe that 
there are large bands of k for which the power spectrum is almost scale invariant. Note that 
there is a feature in the power spectrum on scales k tTa , n that went ultralocal at around the 
time of the transition from expansion to contraction. Also note that modes on larger scales 
(which went ultralocal in the dark energy phase) have a comparable amplitude to those on 
smaller scales (which went ultralocal in the ekpyrotic contraction phase). 

The range in In k of modes to the left of the feature is of order A*d ar k, the number of 
e-folds of dark energy expansion, while the range of modes to the right of the feature is of 
order N ekp = In a/— V^JVo- As we explain later (Sec. |Vj), the value of k for modes on our 
current horizon scale is approximately A e k p /2 e-folds to the left of fc en d and hence near the 
middle of the approximately scale-invariant region of the power spectrum. 

IV. MATCHING PERTURBATIONS THROUGH THE BOUNCE 

Now we map perturbations through the singularity. First we relate physical length scales 
on either side of the bounce. Second, we relate the Newtonian potential for each mode on 
either side of the bounce. 

The first task is simplified by the near symmetry of the contracting and expanding ki- 
netic phases. Consider the wavelength of the last mode to go ultralocal during ekpyrotic 
contraction, with wavenumber fc en d in Fig. 01 Its physical wavelength is roughly given by 
the Hubble radius at that time, ~ 1/V — Knd- By the symmetry of the kinetic phases, its 
physical wavelength when reaches en d on the way out after the bounce will again be 
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FIG. 3: A plot of the power spectrum P$ of quantum-generated fluctuations in the Newtonian 
potential $ going into the crunch. A horizontal line corresponds to scale invariance. /c e nd indicates 
the wavenumber of the last modes to go ultralocal during ekpyrotic contraction. ko indicates the 
scale corresponding to our horizon today. The parameters chosen in Eq. © are b = 0.1 and c = 30. 

For the second task, we need to know how the Newtonian potential behaves on either 
side of the bounce in the four- dimensional effective treatment. All modes become ultralocal 
as the bounce is approached, and so all behave in the same manner, independent of their 
wavelength. If t = corresponds to the bounce, the Newtonian potential goes like A+B t -4 / 3 
on the way in and as A' + B 1 1 -4 / 3 on the way out (see Sec. IVII|) . So each side has a diverging 
term and a constant term. On the way in, the constant term is the "decaying" mode while 
the t -4 / 3 term is the growing mode. On the way out, the roles are reversed with the constant 
term now the "growing" mode and the t -4 / 3 term now the decaying mode. 

The perturbations are in the growing mode approaching the bounce. It is essential for the 
success of the cyclic model that such perturbations lead to some growing-mode perturbations 
after the bounce. For the shape of the spectrum to be preserved it is also essential that such 
mapping occurs in a manner which is at least approximately independent of wavelength. 
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To perform the mode matching one must move from the four dimensional effective theory 
into five dimensions (where the bounce now corresponds to the extra dimension momen- 
tarily contracting to zero size and then expanding again). The work of Tolley, Turok and 
Steinhardt [8] does this and provides us with a matrix -Mtts that relates the incoming 
and outgoing mode coefficients. The form of this matrix is presented in Sec. IV11[ after we 
have introduced normalized mode functions. According to the matching prescription of 
an incoming growing mode maps onto an outgoing solution with a non-zero growing mode 
component. Furthermore, the matching is independent of wavelength, as required. The 
procedure is to find a quantity that behaves as a massless scalar field in a particular five 
dimensional gauge near the bounce. Earlier work of Tolley and Turok 30| showed that there 
is a natural way to analytically continue such fields through the bounce. By understanding 
the correspondence between four and five dimensional perturbations one can thus match 
four dimensional modes across the bounce. Note that alternative matching prescriptions 
with alternative matching matrices are straightforward to use in place of .Mtts- (Indeed, 
very recent work of McFadden, Turok and Steinhardt j3] takes a wider five dimensional 
view of the vicinity of the collision, and effectively pre- and post-multiplies Mtts by another 
matrix. However this does not alter the qualitative features of the matching, so Mtts is 
used in this paper.) 

The result of the matching is that after the bounce the "primordial" power spectrum 
for the expanding phase is in the "growing" mode and has just the same shape as the 
power spectrum before the bounce (i.e. Fig. El again). The amplitude of the power spectrum 
of the Newtonian potential $ is now roughly constant and is determined by three things: 
the starting amplitude, given by the adiabatic vacuum assumption; the amount of growth 
occurring during contraction; and finally the precise matching coefficient of growing mode 
to growing mode in the matching matrix. The parameters controlling the length of the 
contraction phase (e.g. en d) and the details of the bounce (such as the relative speed of 
the branes at collision) must be chosen to make the perturbation amplitude approximately 
10~ 5 , in order to match observations. 
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FIG. 4: A timeline of the cyclic universe showing the behavior of key quantities over the course 
of a cycle. The labels u k tTan exits" and "/c en d exits" show the times where these two modes start 
to follow ultralocal evolution (i.e. when spatial gradients become negligible). 

V. CLASSICAL EVOLUTION OF PERTURBATIONS 

The "primordial" power spectrum computed at the beginning of an expansion cycle can 
be evolved straightforwardly through to the end of the matter epoch in order to compare 
it with observation. In this section, we track the perturbations further around the wheel 
diagram into the dark energy and contraction phases, in order to see what effect they have 
on the quantum generation of the next round of perturbations. We solve the perturbed 
Einstein, fluid and scalar field equations numerically mode by mode, starting deep within 
the radiation era. (All modes of interest follow ultralocal evolution in the kinetic era, so we 
need not worry about the evolution there.) Our numerical code employs synchronous gauge, 
but we express the results in terms of the fully gauge-fixed Newtonian gauge potential 

The perturbations are initially set in their adiabatic growing mode when they are all 
ultralocal (i.e. outside the effective horizon). Observations of the cosmic microwave back- 
ground and large scale structure are sensitive to scales from our cosmological horizon down 
to roughly ten e-folds in k smaller. To relate this to our simulations, we need to know 
what portion of the "primordial" spectrum is relevant to observation. To do this, we work 
out the difference in k between modes on our horizon today and the last modes generated 
during ekpyrotic contraction. As mentioned in the previous section the physical wave- 
length of the latter modes is roughly 1/ \f—V cn( \ when passes en d in the expanding phase. 
Since then there has been a brief w ^> 1 phase, a second kinetic energy dominated phase, 
and the radiation and matter phases (see Fig. Q), providing an expansion of e 27 / 3+7Vrad be- 
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tween them. So these modes have a physical wavelength today of order e 27 ' 3+7Vrad / V — Knd- 
Modes on the horizon today (having a physical wavelength of order 1/vVo) thus have a 
wavelength e~ 27 ' 3_iVrad+iVek p times this, where iV e k p is as introduced previously in Sec. IIII1 
The universe reheats to a moderate temperature after the bang. The reheat temperature 
T rh is tuned to produce density perturbations of the observed amplitude and to satisfy 
other constraints. This imposes a constraint on 7, as defined in Sec. 10] via the relation 
Trh — e_7 ( — Knd) 1//4 (which implies iV rad = — 7 + N ekp /2): namely one needs 7 ~ 10 — 20 (or 
e~ 7 ~ 10~( 4 ~ 8 )), depending on the value of c. For a more precise discussion of the permitted 
range for 7, see Ref. Thus, the modes on the horizon today have a wavelength that 

is e~ 27//3 ~ iVrad a/— V cn d/Vo = e 7//3+iYckp//2 times that of the last modes to be generated in the 
ekpyrotic phase; the former are a factor of 7/3 + > 5 + N ekp /2 lower in Ink than 

the latter. Since iV e kp is very large, of order 100, modes on our horizon today lie roughly in 
the middle of the logarithmic k range of modes to the right of the transition feature seen 
in Fig. |21 The power spectrum is very smooth in the relevant 10 e-folds in k around this 
point, with only a slight tilt. In Table |1J we present a table of some of the scales mentioned 
in this paper and give a timeline of the model in Fig. HJ 

As an example displaying the qualitative behavior of the Newtonian potential, we have 
studied the case for a potential of the form given in with b = 0.1 and c = 5. The code 
starts with seven e-folds of radiation domination remaining. This is followed by seven e- 
folds of matter domination and then only two e-folds worth of dark energy domination before 
ekpyrotic contraction begins. These parameters are not far from those that might give a fully 
realistic description of the universe, and serve to illustrate the qualitative features of the 
mode evolution without requiring us to introduce large exponential factors that complicate 
the numerics. 

Fig. El presents the Newtonian potential at four different times. Note that, in this plot, 
a horizontal line corresponds to scale invariance. Moreover, the tilt of the input power 
spectrum has been neglected, and thus Fig. El depicts the transfer function for an exactly 
scale-invariant power spectrum of the Newtonian potential. 

In curve A, the power spectrum P$ is shown at a time corresponding to "today", the 
end of the matter epoch and the beginning of the dark energy epoch. The spectrum is scale 
invariant at small k and has the expected bend at a scale corresponding to the horizon at 
matter-radiation equality. At larger wavenumber, the curve falls off as k~ 4 . The scalar field 
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TABLE I: A table showing various scales discussed in this paper relative to today's horizon. Note 



that 7 = ln((-Vend) 1/4 /Trh), ^ckp = In V-Knd/^o, and A^ rad = ln(T rh /T ) = -7 + iV ckp /2. 
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fluctuations have had no significant effect on P$ at this stage. This curve is in perfect accord 



with observations. 

Curve B shows the power spectrum evolved through the dark energy phase to turnaround 
at H = 0. On small scales the power has dropped uniformly as perturbations are diluted. 

Curves C and then D show the power spectrum at two stages in the ekpyrotic contraction 
phase. On large scales the spectrum grows rapidly, and on small scales it reddens. The 
unstable scalar field contribution dominates the Newtonian potential in this epoch, as matter 
and radiation are negligible. 

VI. INTERFERENCE OF PERTURBATIONS FROM DIFFERENT CYCLES 

We are now in a position to investigate whether the current perturbations interfere with 
the quantum generation of the new. First, we compare the amplitude of the current per- 
turbations with that of the new ones to be generated, neglecting any effect that preexisting 
perturbations might have on the quantum generation of new perturbations. Then we check 
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FIG. 5: A plot of the power spectrum P$ of the Newtonian potential <3? (in arbitrary units) at 
various points in a cosmic cycle, k = 1 corresponds to the horizon at the start of the simulation 
in the radiation era, and Ink ~ —6.6 corresponds to the horizon at matter-radiation equality. A 
horizontal line corresponds to scale invariance. The curve A indicates the power at the end of the 
matter epoch. Curve B indicates the power at turnaround (H = 0). Curve C indicates the power 
some way into ekpyrotic contraction, and curve D indicates the power further into the ekpyrotic 
phase. 

if this approximation is justified. 

We concentrate on scales that are observationally relevant at the beginning of dark energy 
domination in the next cycle. The mode whose wavelength will equal the current horizon 
radius one cycle from now is today on a scale that is a factor of roughly e~^ dark+ '' Vrad+2 ' y / 3 ) 
times the current horizon radius, the inverse of the total expansion of the universe in the 
interim, corresponding to A^ark + -ZVrad + 27/3 e-folds in k. We saw above that the smallest 
wavelength mode produced in the ekpyrotic phase lay roughly iV ra d + 47/3 e-folds within 
the current horizon. So if there is a significant dark energy phase then there will be negli- 
gible power on the scale of the horizon radius in the next cycle. However, as we shall see, 
an extended period of dark energy domination is not necessary for the cyclic model. To 
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prove the point, we will consider the "worst case scenario," that the number of e-folds of 
dark energy domination is negligible. This would put the horizon radius a cycle from now 
near the high-fc downward turn of the "primordial" power spectrum. Since we have only 
worked to logarithmic accuracy in determining scales, let us conservatively assume that the 
"primordial" spectrum is still scale invariant even on these small scales. Including the red 
tilt will only strengthen the argument below showing that there is no obstruction to cycling 
due to build up of perturbations from cycle to cycle. 

We estimate the amplitude of the current perturbations on the relevant scale as follows. 
The amplitude a mode on our horizon scale has now is roughly 10~ 5 . In the next section we 
show that during ekpyrotic contraction the mode grows like 1/t with t here the time to the 
forthcoming bounce. From our numerical results we see that the tilt of the power spectrum 
on subhorizon scales is reddening as we proceed to the bounce. Asymptotically we expect 
the mode amplitude to drop as k~ A based on the following argument: $ from the matter 
is two powers down from scale-invariant, and this sources a perturbation 5<p a further two 
powers down (see Eq. (fTK|) below, neglecting the time derivative terms). Then when the 
scalar field is again dominant, Eq. (jHJ) tells us that $ should now go like 5(f), four powers 
down. 

Thus with ko corresponding to our horizon, the amplitude of the current perturbations 
on a scale k at a time t before the crunch is roughly 

where to is the time from the start of the ekpyrotic phase to the crunch. So the amplitude 
on our future observer's horizon is roughly 

1Q -5 *£ e -4AT rad (n) 
t V ' 

Now let us estimate the amplitude of the new perturbation mode generated on the same 
scale. The quantum perturbations generated are almost scale invariant, and we exploit this 
by actually calculating the amplitude of a mode generated in the dark energy phase. The 
last of these will have been generated just before t , the onset of ekpyrotic contraction. Just 
as in slow-roll inflation, $ will from Eq. (jSJ) then have an amplitude of roughly (j)/H times 
the amplitude of the fluctuation in 5<p, which is H. So $ is roughly before turnaround, 
which is of order H , the Hubble constant today. During the contraction the new mode also 
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grows like 1/t. So at a time t before the crunch, the newly-generated mode amplitude will 
be of order 

Ho'j. (12) 
Comparing (JTTJ) and (fT2|) we see that the time dependence cancels and we need the inequality 

10 -5 e -4Ar rad < Hq ( 13 ) 

to be satisfied in order for the new perturbations to dominate over the old. 
Rewriting e Niad as T r h/T and H as Tq, we obtain a lower bound on T^: 

T rh > KT 5 / 4 ^, (14) 

in reduced Planck units. Thus the reheating temperature need only be a few hundred GeV. 

As shown in Ref. this is not a difficult condition to satisfy. In that reference, it is 
shown that it is possible to have perturbations with an amplitude of 10~ 5 and satisfy all 
other known constraints for a wide span of T r h above a few hundred GeV ranging up to 
10 10 GeV (or more, depending on the value of c). 

So the lack of growth of modes that enter the horizon during the radiation era, leading 
to a k~ 4 drop in power on small scales, combined with a further k~ A drop during asymp- 
totic scalar-field domination, seems easily sufficient to ensure that new perturbations will 
dominate over the current ones for an observer in the next cycle. We have not even had 
to consider other effects that further suppress small-scale power, such as dark matter free 
streaming, in order to reach this conclusion. 

We now check that current perturbations do not more subtly influence the form of new 
ones by interfering with their quantum generation. Perturbations in the scalar field give the 
main contribution to the Newtonian potential during ekpyrotic contraction, since it is the 
scalar field that is dominating the matter content of the universe during this phase. The 
scalar field perturbations satisfy: 

8<t> + 3H5(t> - = -V H 5cj> + 4$0 - 2$V (j ,. (15) 

CL 

The $ terms on the right hand side of (J 15)) provide the opportunity for current perturbations 
in the Newtonian potential to influence the generation of the new ones. We need to check 
that their contributions are much less than that of the V^Scj) driving term. A quick way 
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to do this is as follows. We shall see in the next section that in a growing ultralocal mode 
$ goes like H/a. We then use Eq. (JSJ) and the background equations to deduce that 5(fi 
behaves like —(f) /a. In the background solution both H and <fi g° hke 1/t, so 8(f) in Planck 
units is roughly equal to $. Hence we need to compare times the newly generated $ 
to Vj, times the pre-existing $. Now and are comparable in Planck units, and we 
have already seen that the newly generated $ is much larger than the pre-existing $ on the 
scales of interest. Thus the influence of current perturbations on the generation of the new 
perturbations is indeed negligible. 

One might be concerned about the effect of nonlinearities in the matter power spectrum 
on small scales in this discussion. We do not think that this is important however, because, 
even if the matter does go nonlinear, this does not change the typical Newtonian potential 
very dramatically. Furthermore, the scalar field does not couple effectively to the matter, 
only via gravity. So both the Newtonian potential and scalar field perturbations can be well 
approximated by linear perturbation theory. 

Because of the large amount of expansion in the radiation era, perturbations in the cyclic 
model do not build up on a given comoving scale, even without much expansion from dark 
energy. Thus, our assumption in Sec. IHIl that there is a long-lasting dark energy phase, is 
not necessary. 

VII. GLOBAL STRUCTURE OF THE CYCLIC UNIVERSE 

Perturbations generated during the dark energy phase start off with an amplitude of 
order the Hubble radius during the dark energy phase, say 10~ 60 . They are amplified during 
ekpyrotic contraction. After passing through the bounce they must have an amplitude of 
order 10~ 5 in order to match observations. Hence there must have been a net amplification 
of the order of 10 55 on the largest scales! Nothing in this argument for the amplification is 
unique to quantum-generated perturbations: classical perturbations are amplified as well. 
Since there are no dynamical effects able to suppress power on scales that never enter 
the horizon, we are forced to conclude that perturbations generated two or more cycles 
ago should today have an amplitude formally far in excess of unity. In this section, we 
investigate these exceedingly long wavelength perturbations in more detail and discuss their 
implications for the global structure of the cyclic universe. 
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In order to track down this large amplification, we first need to understand how a general 
ultralocal Newtonian potential perturbation evolves in time. To start with we sketch a very 
general derivation 3l( of the two linearly independent solutions to the ultralocal perturbation 
equations. This derivation also provides a clue to the interpretation of such perturbations. 

One takes the unperturbed Friedmann- Robertson- Walker (FRW) metric, Eq. and 
considers a (small) coordinate transformation. This changes the metric according to the 
Lie derivative. One then demands that the coordinate transformation is such that the new 
metric takes the Newtonian gauge form, Eq. (|7ji. This restricts the form of the coordinate 
transformation allowed, and it turns out that the most general form of $ so induced is: 

$ = a(i-— [ dt'a(t')) + B-, (16) 

where A and B are constants. Furthermore, the corresponding induced stress-energy tensor 
is adiabatic and has no anisotropic stress. Now, we allow A and B to become slowly- varying 
functions of position. To the extent that second-order spatial gradients can be neglected, 
which is the root of the ultralocal approximation, and that the universe is adiabatic and 
free from aniostropic stress, we now have the general solution for the Newtonian potential 
on long wavelengths. 

Thus every ultralocal perturbation mode in a cycle can be written as «/$ / (/ = 1,2), 
with the basis functions taken from Eq. ()16|) to be: 

$ x = 1 - — f dt'a{t') (17) 
a Jo 

$ 2 = N 2 — . (18) 
a 

In C& 1 we integrate forward to the time t from the big bang at time zero. N 2 is a dimensionful 
normalizing factor, required since $ is dimensionless, and is helpfully chosen to make $ 2 
unity at the time when = en d after the bang. $ 1 is the "growing" mode and $ 2 is the 
decaying mode after the bang. 

Approaching the crunch, it is useful to pick a new linear combination of $ 1 and $ 2 as 
basis functions, namely: 

tt rt CI 

& = 1 + - / dt'a(t') (19) 
a Jt 

$ 2 = N 2 — . (20) 
a 



20 



Here t cr is the time of the crunch, and N2 is a different normalizing factor to N2, now chosen 
to make $ 2 unity at the time when <fi = <p en ^ again, on the way to the crunch. $ x is the 
"decaying" mode and $ 2 is the growing mode going in to the crunch. Note that we have 
used these formulae in deriving the behaviour of the modes in the kinetic phases near to the 
bounce in Sec. IIVI and in getting the approximation 1/t for for the growing mode during 
the ekpryrotic phase in Sec. IVI1 

Our perturbation ajQ 1 may be equivalently rewritten as Knowing how the two 

sets of basis functions are related, the two sets of expansion coefficients are related via the 
matrix equation a = Net , with the matrix N of the form: 



N= , (21) 




where n = — J * cr dt'a(t')/ ' N2 and e = N2/N2. For a typical cycle, n is very large (~ 

amaxWKax/Wend)) ~ V ~ V ^/ V = ) aIld £ is W SIIla11 (~ e-^-iV dajk -2 7 /3) _ 

As discussed earlier in Sec. IIV| it is very helpful to decompose a perturbation into growing 
and decaying modes near the crunch for the purposes of matching it through the bounce. 
We now have and $ 2 ready for this purpose before the bounce. The algebraic form for 
$ x and $ 2 will again serve admirably for giving the growing and decaying modes after the 
bounce, with t now measured from this bounce and N 2 redefined for this passing of e nd- 
We write these new modes for the next cycle as $ next , and our perturbation after passing 
through the bounce will be written as a" 63 '*^^. We can now give the explicit form for the 
mode-matching matrix using the work of Tolley, Turok and Steinhardt, as promised earlier 
in Sec. |lVl With a ncxt = M TTS a, then: 



M TTS = ( 22 ) 




Here m = y 2 eV^ ( f' cnd where v is the (non-relativistic) relative speed of the branes at colli- 
sion. For a typical model m is small but not particularly so. 

This analysis allows us to follow an ultralocal perturbation forward from one cycle to 
the next; if it is described by the coefficients a in the one, it will be described in terms of 
equivalent mode functions by the coefficients o; next in the next, with: 

^ncxt = MttsNo ,, (23) 
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Now, we can find the most positive eigenvalue of the combined matrix .Mtts-^V and thus 
finally extract the the amplification factor per cycle. With typical values for m, n and e 
as indicated above, this eigenvalue turns out to be approximately mn and is indeed large 
because m is so large. We have thus confirmed the heuristic argument given at the start of 
this section for a large amplification of ultralocal perturbations from cycle to cycle. Requiring 
this amplification be enough to take quantum fluctuations to 10~ 5 is one of the conditions 
on T r h considered in Ref . 28[ , as mentioned in Sec. IVI1 

What are we to make of this amplification for classical perturbations? It certainly seems 
that a global view of a universe cycling everywhere with only small perturbations must break 
down after a couple of bounces. On the other hand, we have seen in previous sections that 
as far as physical observers with their cosmological horizons are concerned, the cycling can 
continue indefinitely. We believe the correct interpretation is that, as time passes by, widely 
separated parts of the universe begin to cycle independently of one another, which precludes 
a global FRW picture for the entire universe. Nonetheless, in any given observer's horizon, 
the universe appears to be FRW with perturbations small enough that this region is able to 
continue cycling. 

Perturbations on scales larger than one Hubble horizon still have a small amplitude before 
the onset of ekpyrotic contraction. However, they simply correspond to a small change in 
that observer's background FRW model. So by "recalibrating" the background model before 
ekpyrotic contraction, all superhorizon perturbations can be removed. Power on subhorizon 
scales will be practically unaffected by this change. This "recalibration" might lead to small 
amounts of space curvature in the new FRW background, but this has negligible effect on 
the cyclic history jll|. 

To show that ultralocal perturbations around some point just correspond to a change in 
the background model, we effectively invert our above derivation of the ultralocal behavior of 
the Newtonian potential. Around a chosen point, both the value and first spatial derivatives 
of the real-space Newtonian potential can be set to zero with a dilatative gauge transforma- 
tion. Furthermore, the anisotropic second spatial derivatives can also be removed, leaving 
one as claimed in a different FRW universe with perhaps some modest amount of spatial 
curvature corresponding to the isotropic second spatial derivatives. 

Since different patches suffer different dilatative gauge transformations and then have 
independent Fourier-expanded perturbations, it is clear that we should not expect to be 
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able to sew them back together again at later times and recreate a single global FRW 
solution with small Fourier perturbations; widely-separated parts of the universe are cycling 
independently and out of synch. 

VIII. HOMOGENEITY, ISOTROPY AND FLATNESS WITHOUT INFLATION 

The previous sections have emphasized the importance of the ekpyrotic phase and the 
overall expansion of the universe over the course of a cycle (as illustrated in Fig. EJ) in 
understanding the behavior of perturbations. We now examine their role in explaining why 
the universe is so homogeneous, isotropic and flat today. 

When the cyclic model was originally introduced Q], it was thought that the dark 
energy phase played the critical role in making the universe homogeneous, isotropic and 
flat, as well in ensuring that the cyclic solution was a stable attractor. This being the 
case, some would argue that the cyclic model should rightly be regarded as a variant of 
the standard inflationary scenario, since the dark energy phase can be viewed as a period 
of very low energy inflation. However, as the cyclic model has become better understood, 
we have learned that the dark energy phase plays only a supplementary role in smoothing 
and flattening the universe. In fact, as we shall now explain, homogeneity, isotropy and 
flatness can all be achieved even without the dark energy phase, as was first suggested by 
the "cosmic no-hair theorem" proved in Ref. j^]. 

First, it is already clear from Sees. V-VII that dark energy is not needed to make the 
universe homogeneous. If it were, we would have had to impose the condition that iV dark 1, 
analogous to the condition that the number of e-folds of inflation must satisfy A^ n fl at ion ~ 
50. In actuality, we explicitly assumed iVd ar k = 0(1) and showed that the universe is 
homogeneous after each bang when one takes into account the slowly contracting ekpyrotic 
phase. 

For the curvature, we need to track what happens to fix oc 1/ (aH) 2 . During the ekpyrotic 
phase with w 3> 1, a shrinks by a small amount e~ 2Nckp ^ 3( - 1+w ^ but H grows by a huge factor, 
e -Nekp. go ne |- e ff ec t is that fix is suppressed by a factor of roughly e 2Nckp . During the 
contracting, kinetic energy dominated phase with w ~ 1, and the subsequent expanding 
kinetic and radiation-matter dominated phases, aH undergoes a net shrinkage by a factor 
of e Nvad+41 ^ 3 , as shown in Fig. 4, so £7k is now enhanced by a factor of e 2JVrad+87//3 . From our 
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key result iV c kp ~ 2(iV ra d + 7), though, one finds that the suppression of curvature during 
the contracting phase far exceeds the enhancement during the expanding phase, resulting in 
a net suppression by a factor of e~( 2Arrad+47 / 3 \ a huge net suppression of the curvature even 
for N daxk = 0. This suppression repeats every time the universe goes through an ekpyrotic 
phase. 

For the anisotropy, a similar analysis applies. The anisotropic universe can be described 
by the Kasner metric. The anisotropy in a Kasner universe is characterized by a term in 
the Friedmann equation proportional to a~ 6 . The energy density in 0, though, grows much 
faster during the ekpyrotic phase, as a~ 3 ( 1+w \ If H grows as e^p, then the scale factor a 
shrinks by a factor of e^ 2Nck p^ 3( - 1+w ^ . Hence, the ratio of anisotropy to the scalar field energy 
density shrinks by a net factor of e~ 2iVck p( 3li '~ 1 )/( 3 ( 1+u ')) or roughly e~ 2Nckp in the limit w 1. 
During the kinetic contracting and expanding phases, the ratio is fixed. During the radiation 
and matter dominated phases, the anisotropy is only further suppressed. Hence, like the 
inhomogeneity and curvature, the anisotropy undergoes a net exponential suppression during 
each and every cycle. 

Recall that, in standard big bang cosmology, a puzzling aspect of the large-scale homo- 
geneity and isotropy of the universe is that distant regions within the observable horizon 
were not causally connected in the past. Inflation addresses this aspect by rapidly stretching 
a tiny, causally connected region by a huge exponential factor. In the cyclic model, causal- 
ity is not an issue in the first place because the region that evolved to form the observable 
universe today was only a few meters or kilometers across during the previous cycle, easily 
small enough to have been in causal contact with itself during the previous radiation and 
matter dominated phases. However, in principle, the universe could be causally connected 
and still not be homogeneous, isotropic and flat on large scales. What we have shown in 
this paper is that the ekpyrotic phase and, in particular, the relation iV ck p 2(iV rad + 7) 
automatically insures that it is. 

IX. CONCLUSIONS 

In this paper, we have examined the generation and evolution of perturbations over 
many cycles. First, we have shown that the ekpyrotic phase suffices to make the universe 
smooth, isotropic and flat on large scales. As for the perturbations, we have explained 
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how the galaxies and large scale structure in any given cycle are generated by the quantum 
fluctuations in the preceding cycle without interference from perturbations or structure 
generated in earlier cycles and without interfering with structure generated in later cycles. 
Furthermore, we have examined the global structure of the cyclic universe. Although the 
universe can be described as a nearly uniform Friedmann-Robertson- Walker within any 
observer's horizon, we find that global structure is more complex and cannot be characterized 
by a uniform Friedmann-Robertson- Walker picture. 

An important corollary of our results is that neither an extended dark energy phase nor 
any other form of inflation is needed to solve the horizon and flatness problems. Instead, 
the universe is made sufficiently smooth, isotropic and flat during each ekpyrotic phase in 
which the universe contracts with w ^> 1. Hence, the cyclic model should not be construed 
as a variant of inflation. Rather, the ekpyrotic contraction mechanism should be viewed as 
a genuinely novel approach for solving the classic cosmological problems. 

A fuller understanding of the bounce and a fundamental derivation of V{4>) remain the 
most pressing issues for the cyclic model. Major surprises there aside, our work has shown 
that the cyclic model is in good shape as a candidate for a complete cosmology for the 
universe. 



APPENDIX: THE w > 1 EXPANDING PHASE 



In the cyclic model, the scalar field has to acquire a boost to its kinetic energy at or 
after every brane collision, in order to overcome the additional Hubble damping due to the 
radiation and to make it back onto the potential plateau. 

As discussed in the original papers the boost may be parameterized as 



dt 



dt 



out 



(24) 



with x a small parameter. Such a boost can be produced either by the production of extra 
radiation on the negative tension brane, or by the nonminimal coupling of <j) to matter, 
which drives it positive in the expanding phase. 

Both on the way in to the bounce on the way out from it, the energy density of the 
universe is dominated by the kinetic energy of the scalar field. For small \i the outgoing 
solution is nearly the time reverse of the incoming one and, as the field <fi crosses </> e nd ; the 
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solution is close to the time reverse of the scaling solution, Eq. |S| However, the scaling 
solution is not an attractor in the expanding phase, so small deviations from it grow with 
time. The modest increase in scalar field kinetic energy, parameterized by \i causes the 
kinetic energy of to eventually overwhelm the potential energy, so that the solution enters 
a second kinetic energy dominated phase. 

We compute the value of where this second expanding kinetic phase begins by elimi- 
nating t in favor of <j) in the background equations © and Q, obtaining 



dH 3H 2 - V 

^=-v^^- (25) 

Approximating V ~ — Voe~ c< ^, we then change variables to h = He c ^l 2 to remove the leading 
dependence in the scaling solution, getting 



* = £ k - J 3 J^, (26) 
d6 2 V 2 K ' 



from which the fixed point scaling solution h sc = y Vo/(|c 2 — 3) is recovered. Now, we can 
describe the effect of the small perturbation x by linearizing (|2T)j) about the scaling solution, 
obtaining 

d5h fc --^]5h. (27) 



2 C/ 

Thus, the small perturbation 5h grows exponentially with 0. The initial conditions for Sh 
are found from the scaling solution, H 2 = <f) 2 /c 2 , the Friedmann equation and the definition 
of h to be 5h/h sc = (c 2 /6)(<50/0) = c 2 x/6. The perturbation grows until Sh/h sc is of order 
unity, when 

-^ d " c(l-6c- 2 ) ln ^' (28) 
after which the potential V((j)) becomes irrelevant, and, from (|25j) or (|26jh H oc e^V^/ 2 "^ 
the expanding kinetic energy dominated solution. If c is large and \ is not extremely small, 
the second kinetic phase starts rather soon after <fi passes e nd- It follows that the w ^> 1 
expanding phase is brief and can for most purposes be safely ignored. 
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